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Nomenclature 


fll ••••« 
l e 

aging time shift factor, at aging 

time t e 


....aging time shift factor, at aging 


time t+t°, relative to aging time t° 

B tenf " 

shift rate at reference aging time, 

unity by definition 

«r 

time temperature shift factor 

-t'j 

fl li/Tj • 

time-temperature shift factor 

between Tj and T 2 at aging time 

[A] 

....matrix relating midplane strains to 
resultant forces in laminate 

P 

....shape parameter 

C°i 

....free volume parameter for WLF 
equation 

C°2 

....free volume parameter for WLF 
equation 

{ £ *} •• 

....vector of strain components = 

{f xx ' e yy' e xy } 

Kl- 

....strains at the midplane of the 
laminate 

h 

....thickness of laminate 


....vector of curvatures at the 
midplane of the laminate 

X 

....effective time 

M 

....shift rate 

Ml 

....shift rate prior to equilibrium 
aging time 

M2 

....shift rate after equilibrium aging 
time is reached 

M 22 ••••■ 

....shift rate for S 2 2 

M66 

...shift rate for S 66 

{M,}. 

...resultant moments in laminate 

{N,}- 

...resultant forces in laminate 


[P(0].... effective modulus matrix of the 

laminate in the loading coordinate 
system 

0 angle between the fiber direction 

and the longitudinal loading axis 

[01... .. inverse of the S-bar matrix for the 
kth lamina 

S creep compliance 

So initial compliance 

[S ] effective compliance matrix 


S 22 

... transverse compliance in 0° lamina 

S 66 

...shear compliance in 0° lamina 

<*>. 

...constant stress in time 

i a x }•“ 

... vector of stress components = 
fyxx’° w a xy } 

fix 

... average stresses through the 
thickness in laminate 

t 

...time 

T 

... relaxation time 

te 

... aging time 

«<” 

... initial aging time 

*CT€f *••• 

... reference aging time 

* 

*c 

...equilibrium aging time 

m 

... transformation matrix based on 6 

T i 

...glass transition temperature 

Tref- 

...reference temperature, where C°, 
C® are defined 

2 

... coordinate through the thickness of 
the laminate 
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Abstra ct 

For many polymeric materials in use below the glass transition temperature, the long term 
viscoelastic behavior is greatly affected by physical aging. To use polymer matrix composites as 
critical structural components in existing and novel technological applications, this long term 
behavior of the material system must be understood. Towards that end, this study applied the 
concepts governing the mechanics of physical aging in a consistent manner to the study of 
laminated composite systems. Even in fiber-dominated lay-ups the effects of physical aging are 
found to be important in the long-term behavior of the composite. 

This paper first lays out in a self-consistent manner the basic concepts describing physical 
aging of polymers. Several aspects of physical aging which have not been previously 
documented are also explored in this study, namely the effects of aging into equilibrium and a 
relationship to the time-temperature shift factor. The physical aging theory is then extended to 
develop the long-term compliance/modulus of a single lamina with varying fiber orientation. 
The latter is then built into classical lamination theory to predict long-time response of general 
laminated composites. Comparison to experimental data is excellent. In the investigation of fiber 
oriented lamina and laminates, it is illustrated that the long term response can be counter- 
intuitive, stressing the need for consistent modeling efforts to make long term predictions of 
laminates to be used in structural situations. 


Thomas S. Gates 
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MS 188E 
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Introduction 


The High Speed Civil Transport aircraft (HSCT), envisioned to have a lifetime of over 60,000 
flight hours and to travel at speeds in excess of Mach 2, is the source of intensive design and 
development studies at NASA and major airframe developers. In the area of materials for the 
HSCT, polymer matrix composites (PMC’s) are being strongly considered for use in primary and 
secondary structures due to their high strength to weight ratio and the options of property 
tailoring. However, one difficulty in the long term use of polymer based materials is that their 
strength and stiffness properties may change measurably over time (especially at the elevated 
temperatures experienced during flight). Analytical prediction of these properties based on 
irregular thermal and mechanical loading may be extremely difficult. 

The current study is focused on one aspect of long-term polymer composite viscoelastic 
behavior: physical aging. In contrast to the equilibrium state of a polymer above its glass 
transition temperature, when a polymer is cooled to below its glass transition temperature the 
material does not achieve instantaneous thermodynamic equilibrium. Instead the free volume, 
enthalpy and entropy evolve over time toward their hypothetical equilibrium values. During this 
time the mechanical properties change significantly. This change is termed physical aging. Figure 
1 illustrates the change in the volume in a polymer as a function of temperature; when the 
polymer is cooled below the glass transition temperature, the volume is initially in a non- 
equilibrium state and evolves slowly towards equilibrium volume with time. 

In essence, during physical aging the material becomes suffer, so that at a given time the 
compliance is decreased (or the modulus increased) over what one would have expected in a 
viscoelastic material without aging. In terms of free volume theory (Struik, 1978; Ferry, 1980), 
one can visualize that as the free volume decreases towards its equilibrium value, the mobility of 
chain segments is hindered, giving rise to a stiffer response. Aging is a characteristic of the glassy 
state and is found in all polymer glasses, as thoroughly documented by Struik (1978). The aging 



of polymer composite systems has received little attention until recent years. But several 
thorough experimental studies (Sullivan, 1990; Gates and Feldman, 1993; Hastie and Morris, 
1993; Sullivan, Blais et al., 1993) have illustrated that the matrix dominated composite properties 
of continuous fiber reinforced PMC’s, namely the shear and transverse response, are affected by 
physical aging in a similar manner to pure polymers. 

This work begins with a review of the concepts of physical aging on a pure polymer system 
and the theory is outlined in a consistent manner. The effective time theory (Struik, 1978; Knauss 
and Emri, 1981), which can be used to predict long term behavior based on short term data, is 
also mathematically formalized. The effects of aging to equilibrium are proven and discussed and 
a special relationship of aging to time-temperature superposition is demonstrated. The physical 
aging theory developed for polymers is then applied first to a unidirectional composite, then to a 
general laminate. It will be shown that the effects of aging on the long-term properties of 
composites can be counter-intuitive, stressing the importance of the development and use of a 
predictive theory to analyze PMC structures. 

Physical Aging of Polymers 

The traditional means to assess and quantify the effects of physical aging on mechanical 
properties is through a series of sequenced creep (constant load) and recovery tests (Struik, 
1978). In these tests, the specimen is initially quenched from above the glass transition (T g) to a 
temperature below Tg. The time the material exists below its glass transition is referred to as the 
aging time, t e . As aging time progresses, a series of short (in comparison to the elapsed aging 
time) creep tests are run to extract the momentary creep compliance of the material. The 
momentary creep compliance of the material can be described with a three parameter fit model, 
first proposed for linear viscoelastic polymers by Kohlrausch (1847), 

Sit) = S(je (t t r)P (1) 

where S 0 , is the initial compliance, j9 is a shape parameter, t is time, and x is the relaxation time. 



As shown in Figure 2a (and more explicitly in the isochronal plot in Figure 2b), the material 

is stiffer (or less compliant) at any given test time as aging time progresses. In nearly all cases of 

aging of polymers and PMC’s it is possible to bring the momentary creep curves measured at 

different aging times into superposition through a horizontal shift.* This result is akin to the well- 

known time-temperature superposition principle (Ferry, 1980). If the aging time shift factor, a t , 

is plotted as a function of aging time on a double-log scale, it is found to map a straight line with 

a slope of /J. as shown in Figure 3. The shift rate, fl, 

dloga t 


V = ~ 


dlogt e 


( 2 ) 


usually has values on the order of unity and can be considered to be a material constant. Thus, a 
series of creep tests (as in Figure 2a) can be used to experimentally determine the value of the 
shift rate, JJ., for any given material. 


Since the momentary creep curves superpose through horizontal shifting on the log scale, 
the only parameter which changes as a function of aging time is the relaxation time. In terms of 
the mathematical expression of compliance in Equation (1), the momentary material properties 
then vary with aging time according to 

S(t>- Soe"^ ( 3 ) 

where 

T(t e ) = T(t eref )/a te (4) 

is the reference aging time, a t =1 and the aging time shift factor is given by 


t A small vertical shift is also sometimes needed in order to achieve ideal superposition. The full implications of 
the vertical shift factor in terms of modeling will be reserved for a later investigation. A study of the vertical 
shifting involved in aging of composites can be found in (Sullivan and Blais et al., 1993). As with all viscoelastic 
material descriptions and shifting procedures, certain trends have been shown to be valid for a wide range of 
amorphous polymers (for physical aging, Struik (1978) has compiled exhaustive data) and the equations to 
describe the behavior are selected based upon this body of empirical data; the constants used in the equations are 
material properties determined from specific experiments on each polymer individually. 



( 5 ) 


Thus to describe the momentary creep compliance at any aging time, the parameters needed are 
the initial compliance ( S 0 ), the shape parameter (/?), the shift rate C/i), and the relaxation time ( z) 
at a reference aging time ( t ere f). Note that if a given aging time is greater than the reference aging 
time (t e > t ere f), a t( < 1; otherwise a te > 1. In the case that t e >t ere f, the relaxation time at t e is 
also greater (r (t e ) > r (t ere f)) and for any given test time the material is in a “glassier” state such 
that the compliance is lower and the modulus higher ( S(t; t e ) < S( t; teref»- 
Effective Time Theory 


If the loading time of a material is not short in comparison to the previous aging time, then 
aging proceeds even as the test progresses. Thus, if the beginning aging time is f® , at some time 
later during the test the aging time is f + ff, where t is the test time. When the test time, t, 
approaches and exceeds the value of the initial aging time, one would expect deviation of the 
long-term creep response from the momentary. Taking the initial aging time t° e to be the 
reference aging time ( t e ref = tf), the shift factor, a t o , at any instant in time can be defined based 


on fi -. 


a t o(t) = 


( t o Y 

l e 

t°e+t) ‘ 


( 6 ) 


In this notation, the updated aging time is always greater than the initial aging time and 
consequently a t o <1. 


To account for the cumulative effects of aging, one can consider what happens in any 

discrete differential time interval, dt. During a time increment dt at time t all relaxation times are 

1/a slower and the time increment is 1/a less eventful than if aging time had not been increasing. 

The “effective” time increment can then be defined (Morland and Lee, I960; Struik, 1978) 

dX = UfO (t)dt 

and the total test time can be reduced to the “effective time” (or “reduced” time) , A 
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( 8 ) 


X = fa,.(4>dl 

0 


Integration of (8) using (6) leads to two distinct expressions for effective time 

li = l ... X = t° e ln(± + 1) 


... A = 


i-n L 


(2 + 




(9) 


Using the effective time in place of real time in Equation (3) 

»D=V' 1miV «„) 

results in prediction of long-term material response (as opposed to the momentary response) 
based only on viscoelastic material parameters determined from short term tests. Notice that the 
relaxation time, r, is calculated based on a reference aging time, t^, which may differ from t° e , 
in which case the value of t is updated on the basis of equations (4) and (5) using t e mt° e , not 
equation (6): 


f 

T(t?) = T(t eref ) 

V 



^eref J 


(ID 


At any given test time, f, the long term response is mapped from the value of the 
momentary response at time A. Figure 4 illustrates the effective time theory prediction of long- 
term material response using idealized material properties. Note that the compliance deviates 
sharply from the exponential momentary curve as the test time approaches and exceeds the 
initial aging time.’ If the initial aging time for a long-term test is changed, both A and r must be 
updated based on the shift rate and the new initial aging time, using equations (6) and (9) and 


For the illustrated long-term test, the material response would follow the long-term response curve only, not the 
momentary curve. Experimental data can be obtained for the upper portion of the momentary curve, however, by 
performing tests at higher temperatures with the same initial aging time and using time-temperature superposition 
to get the momentary master curve. (At.the higher temperature, the upper portion of the momentary curve is 
actually obtained at shorter times, ie, before f > tg.) 
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(11). For n=I, the long-term curves for different t° are simply shifted horizontally from one 
another on the log time scale; for fi*l this is no longer true. 

Figure 5 illustrates the role of the shift rate, fi, in the long-term response of a material. Note 
that increasing values of fi lead to a stiffer (less compliant) response at longer times and that a 
small change in the shift rate can lead to radically different response at long times. Thus the shift 
rate can be used as a viscoelastic screening parameter for selection of materials. As the shift rate 
goes to zero, the long-term response approaches the momentary response, so that a shift rate of 
zero indicates no aging and the long-term curve coincides with the momentary curve for t ^ ; the 
creep compliance then increases exponentially with loading time without the retarding effect of 
aging, a situation undesirable for most PMC structures. As a final note on the significance of the 
shift rate, in the past shift rates were considered to have an ideal value of unity, but values of less 
than one were also common. No mention was made of shift rates exceeding one, particularly in 
the theory, but current data on composites has yielded shift rates exceeding one (Gates and 
Feldman, 1993). In this respect, it should be noted that the aging and effective time theory is self- 
consistent and shift rates exceeding one are both physically and mathematically permissible. As 
the shift rate increases from a value of unity, the effective time A decreases and if the shift rate 
were allowed to obtain the physically unrealistic value of infinity, the effective time would be 
zero and the compliance of the material would equal the initial compliance at all times. 

Aainq to E quilibrium 

For most polymers, if the material is being loaded at temperatures close to its glass 
transition, the material can age into effective equilibrium during a relatively short time frame. The 
time required to achieve equilibrium is known as the equilibrium aging time, t e (Lee and 
McKenna, 1990; Santore, Duran et al., 1991). Because the material state upon cooling deviates 
further from thermodynamic equilibrium as the temperature decreases and because aging 
proceeds on a log-time basis, the temperature must be very close to T g for equilibrium to be 



reached in a realistic time scale. Such situations are of interest for polymer composite design 
studies because viscoelastic material properties of an aircraft wing, for example, may vary 
significantly if equilibrium is reached. 

According to experimental evidence on a polymer by (Lee and McKenna, 1990) and 
preliminary evidence by (Gates and Feldman, 1993) on a composite, when the equilibrium aging 
time is exceeded the shift rate does not become zero but decreases by an order of magnitude. 
This might indicate that aging has not completely ceased, such that absolute thermodynamic 
equilibrium has not been attained. But recalling the sensitivity of the long term modulus to the 
shift rate, an order of magnitude change in n indicates that the rate of aging is indeed greatly 
reduced. In general, the experimental results showed /Xj « / before £ and n 2 = 0 J after £ ancl 
this is illustrated by the bilinear shift rate in Figure 6. 


It will be shown that the basic concepts of aging and effective time can be mathematically 
formalized to handle the case of aging into equilibrium. The analysis must simply proceed along 
lines considering the distinction between an initial aging time that is less than versus one that is 
greater than the equilibrium aging time. To calculate the momentary compliance at any aging 
time, for t e < £ all precepts remains as before and 


V*i 


T(t e ) = T(t eref ) 


(W J 


( 12 ) 


If t e >£, however, then the calculation of the shift factor (and therefore the relaxation time) must 
account for the change in shift rate after aging into equilibrium. Expressing /x 2 ancl Ml as the 
slopes of the double-log shift factor -aging-time curve (Figure 6), 


lOg\ ~lOg a t; 1 °%K ~ l 0 Z a t'rrf 

log t e - log t\ 1 logt] - logt eref 


(13) 


and combining and simplifying these equations (noting that log a t = 0 by definition), one finds 


l °g\ 


/y v*r 

2 . 

J 


= loga te -logfe-l 


V 


(14) 


S- 



From equation fl 4) the shift factor for L > L is obtained 


a t, = 


( *Y2 

ft ^ 

1 eref 





* 

^ j 


( 15 ) 


from which one deduces that relaxation time can be expressed as 


x(t e ) = x(t ere f) 


( L Y* 


\ l e J 


( .* Y' 


*eref J 


( 16 ) 


In the calculation of effective time for long-term response predictions, the situation is somewhat 
reversed. If t e < t e , then as the test time increases in a long-term test, the aging time increases 
from its initial value of t° e and at some point will exceed the value of t * . Define this point in time 
as t 

t = ~ . (17) 

The integration for the effective time calculation from equation (8) must be separated then into 
two parts. For t < t , the shift factor is as given for equation (6), 


( t o Y*i 

l e 


a t (t<t) = 

l e 

but for t > t the calculation must proceed along the lines as for equation (15), yielding 

a t (t>t) = 


(18) 


t 1 
l e 


f,o\ 

[e_ 

t° + t 
\ l e + l J 


* 

) 


(19) 


Combining, the effective time for i e < t e is 

.o V* * 

where a t o(tg) = l. This can be solved explicitly to yield 


“H it, 


l e 


te 




( 20 ) 


v V ) 


r 

A = — 


1 + k — k 

.0 

\ te J 


-1 


( fO Y* 
\t e ) 


1-V 2 


* 

V te J 


-1 


( 21 ) 


If t°>t e , the effective time for long-term response can be calculated much more readily. 
Again defining a t o (t°) = l, Equation (8) applies directly and 


- 9 - 




( 22 ) 


Equations (12) and (16) can be used quite accurately to predict the momentary response of 
the aging data of a polymer or a composite aged to equilibrium (illustrated in Figure 2a, where 
the curves coincide with experimental data points). Figure 7 shows the results of the long-term 
predictions using equations (21) and (22) along with the momentary data from (Lee and 

McKenna, 1990). 

In Figure 7, note that if the material is loaded with an initial aging time less than t e , the 
response initially follows the momentary curve; then begins to deviate (becoming much stiffer 
than the momentary curve predicts) as the test time approaches initial aging time; then as the test 
time and initial aging time combined exceed t ’, the material has “aged into equilibrium” and the 
response follows parallel to the momentary curve once again. Since the material is aging at a 
much reduced rate after t * , accumulating time no longer has the effect of causing the material to 
remain much stiffer than the momentary response. In the long term response of a material aged 
into equilibrium (t° > (*), the compliance increases in an exponential manner. The results are in 
good agreement with the short term responses at increasing aging times, however there is as yet 
no experimental data for long term response. From short-term data as shown in Figure 2a, after 
equilibrium the momentary responses do not continue to shift significantly on the horizontal log- 
time scale. Thus the long term response, which is in essence a compilation of momentary 
responses at longer and longer aging times, would also not be expected to have the compliance 
values extend horizontally at a consistent low magnitude across the time scale. These results 
indicate that due to the dramatic change in material response upon aging into equilibrium, this 
possibility should be considered during the design stage for structural components of polymers 
or PMC’s, particularly in high-temperature applications. 
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Aging Ti me - Temperature Relationships 


In the course of this study the relationship between aging time and time-temperature 
superposition was investigated. In general, to obtain the most accurate results for the aging time 
shift rate it is necessary to run momentary creep tests at various temperatures for each aging time. 
In this manner, through time-temperature superposition one obtains a larger segment of the 
momentary master curve and the comparisons of momentary master curves at different aging 
times to calculate the aging time shift rate is more accurate. A general time/ temperature/ aging- 
time shift factor can be expressed (Sullivan, 1990) 

loga = loga te +loga T> (23) 


where is the aging-time shift factor as defined by equation (5) which brings isothermal creep 
curves at different aging times into superposition and a T is the time-temperature shift factor 
which is generally expressed through the (Williams-Landel-Ferry) WLF equation at temperatures 
above the Tg (Ferry, 1980) 


__ cfcr- y 
° gaT cf+r-7^ 


(24) 


where Cj and Cf are material constants based upon free volume. In (24) the bar is deliberately 
omitted from the a T since above the T g the time-temperature shift factor will be independent of 
physical aging, as aging does not occur. Unfortunately a theoretically derived equation to 
represent a T below the glass transition temperature is not available, with or without 
consideration of physical aging effects. A common technique when using time-temperature 
shifting on experimental data both above and below the Tg is to use a WLF fit for data above the 
Tg and to do a linear fit for the data below the Tg (Rusch, 1968; Brinson, 1990). 


It has been demonstrated that the shift rate which characterizes the aging of a material, 
although constant for a wide range of temperatures, is in fact a function of temperature (Struik, 
1978). In particular near the Tg of a material, the shift rate decreases, the value being identically 
zero above Tg. Because of this temperature dependence of /r, and thus of a tg , then that the time- 
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temperature shift factor, a T , must be a function of aging time below Tg. It will be shown here that 
exactly this relationship will allow us to mathematically formalize the dependence of a T on aging 
time, allowing us to predict a T at any aging time if we know its value or functionality at a single 
aging time and we know the temperature dependence of fi. 


To see the relationship between a T and fl(T), consider the following four cases outlined 
schematically in Figure 8: momentary compliance curves at temperature Tj for various aging 
times, momentary compliance curves at temperature T 2 for various aging times, momentary 
compliance master curves at t e i for various temperatures, and momentary compliance master 
curves at t e i for various temperatures. Taking t el =t ere f, at Tj the relationship between 
compliance at two aging times can be expressed (using equation 5) 


Sj'(t) = sj‘ 

l el l e 


Similarly at T 2 , 




f ( t \^\ 

ill 

Jelj 


{ t 


(25) 


e2 


lel 

\*elj 


(26) 


For constant aging time t e j, the relationship between momentary compliance master curves at 
two temperatures can be given by 

Kb) (27> 

Similarly at tel . 

S f?2 (t) = S l\ ( ta T% ) (28) 

denoting as the time-temperature shift factor between Tj and Ti at aging time t e j. Equating 
like compliance terms using straight-forward substitution and reduction, one can obtain 

uJ 


a l '2 yi( T 2 )-ii(T t ) 

OhSlL J 
•A. 


e ,n 2 


(29) 


Thus, given the shift rate as a function of temperature, fiCO, and the values or an 
expression for the time-temperature shift factor at a single aging time, a T can be calculated at any 
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aging time. Since time-temperature superposition is generally used to collect sufficiently detailed 
momentary curves at any given aging time, these formal relationships between aging time and 
time-temperature superposition will be of utility in data reduction. 

Aqina and Long Term Response of a Composite Lamina (Orthotropic Plate) 


The response of a composite lamina under plane stress conditions is generally described 
through transformation expressions which relate the properties of the lamina measured along 
and transverse to the fiber direction with the angle, 9, between the fiber direction and the 
longitudinal loading axis. (See Figure 9 ) In this manner an effective compliance matrix, [s], for 
any off-axis lamina is defined (Jones, 1975; Halpin, 1992) as 

[sHir'isp] 


where 


m= 


cos 2 9 

sin 2 9 

2sin9cos9 


$n 

$ 12 

0 ' 

sin 2 9 

cos 2 9 

-2sin9cos9 

[S] = 

$12 

$22 

0 

-sin 9 cos 9 

sin9cos9 

cos 2 9 - sin 2 9 


_ 0 

0 

$66. 


(30) 


The terms in the S matrix are the properties of the lamina with respect to the fiber coordinate 
system and relate the stress to the strain in the fiber direction. The elastic stress-strain relation for 
any fiber orientation under in-plane loading is then given by 


(3D 


f "1 

£ XX 


a xx 

,£ yy 

-P} 

a yy ' 

£ xy 


a *y. 


where £,. and cr„ represent the strain and stress in the axis directions indicated. 


Experimental data on polymer composites (Sullivan, 1990; Hastie and Morris, 1993) has 
shown that the transverse compliance, S 2 2 , and the shear compliance, S 66 , are the only terms in 
the compliance matrix that are dominated by polymer matrix properties. Only these terms are 
thus time-dependent and subject to physical aging. The transverse and shear compliance depend 
on the viscoelastic parameters introduced earlier and can be represented by equations (3), (4) 
and (5), such that in functional form 
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$22 (t)—f( $22 > $22 ' ^22 ( teref)> V-22' D 

$66 (*) = f( $66 ' $66 ' T 66 ( K ref )> 1*66 >*K 


( 32 ) 


Note that the shear and transverse compliance of a composite are each described 
independently by the four viscoelastic parameters: initial compliance, shape parameter, a 
relaxation time at a given reference aging time*, and a shift rate; these parameters are measured 
in the laboratory for a given composite material in short-time tests. Using equations (3) or (10), 
(4) and (5) with equations (32) enables the calculation of the momentary or the long-term shear 
and transverse compliance of a composite lamina with respect to loading in the fiber direction. 


In this analysis, we will not consider the full hereditary integral constitutive expression for 
the viscoelastic behavior (Findley, Lai et al., 1989), but will concentrate on the case of creep 
compliance from constant load. A general linear viscoelastic constitutive law may be written as 


£(t)= f S(f — x)-—j——dx 
J dx 


(33) 


and in a constant load creep test the stress may be represented as 

a(t) = o 0 H(t) 

where H(t) is the unit step function. Substitution of (34) into (33) yields 

t 

e(t) = a 0 js(t- r )S( x)dx 


(34) 


(35) 


e(t) = cr 0 S(t) 

Given a constant load assumption, equation (3D can then be applied directly, substituting the 
appropriate time-dependent quantities. 




Most generally, the reference aging time for the relaxation time can differ for the transverse and the shear 
compliance. One must then, however, carry a superscript along on the reference aging for clarity: tj 2 ^ and 
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Thus, the response of a lamina for any angle, 0, of fiber orientation, can be determined 
through the transformations accounted for in equation (31) together with the expressions for S 2 2 
and $66 and the constant terms of the S-matrix: 

; S ij (t)=f(0,S ij ) (36) 

or in shorthand notation: 

{e x (f)} = [S(f)]{<r. t }. (37) 

For linear viscoelastic behavior, the relaxation times and effective time are updated 
individually for each of the time-dependent terms and then combined in a consistent manner to 
yield the overall long-term response. Note again that equation (37) holds only for the case of 
constant loading; however, the time-dependent compliance matrix itself (but not the stress-strain 
relation) calculated in this manner holds for any loading. 

Note that the S-matrix will generally be fully populated and all of the terms will be time- 
dependent, in contrast the S-matrix. Consequently, due to the complex interaction of all these 
terms and the independent time parameters (fl and T) of S 22 and S 66 , the long-term compliance 
of a lamina with a given fiber angle cannot be intuitively predicted. 

To illustrate the complexity of long term lamina response, consider the following example 
using IM7-8320 data (Table 1) (Hastie and Morris, 1993). In Figure 10 the long term load-direction 
compliance of a composite lamina at various fiber angles is presented as predicted by the theory 
outlined in this paper. The results coincide precisely with experimental data (not shown) up to 
106 seconds, after which there are no further experimental data points. Note that at all times the 
90° lamina is the most compliant, the 0° the least compliant, and the remaining angles falling in 
order in between, as one might expect. To easily compare the relative change over time of the 
compliance with respect to initial value for the different angles. Figure 11 shows the same results 
shifted down uniformly to the 0° curve. In this figure, one sees that the long time response does 
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not consistently follow a “simple” rule that the 0° compliance changes the least (none) with time, 
all others increasingly more so to a maximum with the 90 c lamina. 

Tafab 1 : Transverse and Shear Complianoe Data on M7-8320 system (Hasfie and Mortis, 1 993) 


VE parameters: 

^22 

$66 


0.77 

0.93 

P 

0.416 

0.456 

T 

1.19x 10 6 sec 

4.31x 10 5 sec 

So 

750 xlO' 9 Pa" 1 

1364 x 10' 9 Pa -1 

^ eref 

3.24 x 10 4 sec 

3.24 x 10 4 sec 

Elastic parameters: 

S n = 40.23 x 10~ 9 y psl 

v 12 = 0.348 


If one notes that in Figures 10 and 11, the shift rates for S 22 and S 66 are quite different 
([Iqq = 0.93,^22 =0.77'), an immediate conclusion is that since the shift rates are different for 
those two terms, the aging rates and thus effective time aiso differ. Such a discrepancy will be 
manifested in different rates of change of the long-term compliance for different fiber angles, 
where the relative influence of S 22 or S 66 is determined through the transformation equation 
(30). The complex coupling of the extension and shear behavior is clearly illustrated by 
artificially imposing a constraint of equal shift rates. With such a constraint, it might be expected 
that a “simple” rule for change of compliance over time would be met. However, as is illustrated 
in Figures 12 and 13, if Use = /*22 = 0-77 the result is that even if we assume that the rates of aging 
are the same, the long term response at various fiber angles definitely do not follow any simple 
rule. In fact, in this case for very long times, the compliance of a 43° lamina exceeds that of a 90° 
lamina. And if the curves are again shifted uniformly to the 0° curve it is clear that the 
compliance of the 30° lamina in fact changes the most over time and that the 90° changes the 
least of all except for the 0° lamina. 

In this second case (Figures 15 and 16) where we have artificially forced equal shift rates 
(and equal effective time) for the transverse and shear compliance, one clearly sees the influence 
of another time factor. The relaxation times for the S 22 and S 66 compliance terms also directly 
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affect the long-time response of the lamina. Since these relaxation times in general differ between 
S 2 2 and S 66 , and the relative influence of S 22 versus S 66 in the load direction compliance is 
different for every fiber angle via equations (30) and (37), even if the shift rates are the same it 
cannot be expected that the time-dependent response follow a simple pattern. In this particular 
case, the S 66 relaxation time is an order of magnitude smaller than for S 22 at the same reference 
aging time, indicating that the S 66 exponential will have a larger magnitude for lower values of 
time than S 22 - Thus, for hypothetically equal shift rates, a lamina with fiber angle of 30°, which is 
strongly influenced by S 66 , is predicted to change more over time than the 90° lamina which is 
influenced only by S 22 . 

Conversely, when the shift rates are allowed to assume their real values, since /i 66 has a 
higher magnitude than fi 22 the effect of physical aging will be such that the rate of increase of the 
S 66 compliance with time is substantially more dampened than the S 22 compliance. The 
opposing effects of the shorter relaxation time and the higher shift rate then combine to yield a 
response that is closer over the long term for S 22 and S 66 than expected. An illustration of this is 
given in Figures 14 and 13, where the individual compliance terms are plotted for the two shift 
rate cases. 

Aging of a Laminate 

The next logical step in the analysis is to examine the response of composite laminates to 
physical aging. This can be accomplished by incorporating the process of physical aging of a 
lamina (as developed in the previous section) into standard lamination theory. In this study, the 
theory for thin laminates was used and the results presented are for symmetric laminates with in- 
plane loading only, although the theory applies readily to non-symmetric laminates and more 
general loading. From classical plate theory, the stress-strain relation for the kth lamina in a 
composite is given by (Jones, 1975; Halpin, 1992) 

{<*x} k = [Q]*{{ e *°} + z fcj}, (38) 
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where {<J x } k represents the stresses in the kth lamina, [Q^ is the inverse of the S-bar matrix (as 
defined by equations (30) and (37) for the kth lamina, jej} are the strains at the midplane of the 
laminate, z is the coordinate through the thickness and {k* x } is the vector of curvatures at the 
midplane of the laminate. Upon introduction of the resultant forces (N) and moments (M) in a 
laminate containing n lamina with an overall thickness, h, 

h / 2 n 

{ N x} = |K}<fe = Xj* K} k dz (39) 

-h/2 k=l Zl_1 



(40) 


and substitution of equation (38), one obtains a system of equations which can 
determine the state of stress and strain at any point in the laminate: 
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be solved to 


(41) 


For the remainder of this section, we will be concerned primarily with the A-matrix, the 
upper left quadrant in (41), 

h] = Z[Q]/Zit-W, (42) 

k=l 

since in the case of symmetric laminates and in-plane loading, only this portion of the square 
matrix in equation (41) is both non-zero and relevant. In this case, the A-matrix leads directly to 
the calculation of the effective moduli of the laminate. Because we are examining the effects of 
aging without recourse to a general nonlinear viscoelastic constitutive law, the restriction of 
constant stress applied to the laminate still applies. Denoting with an overbar the average 
through the thickness of the stress applied to the laminate, equation (39) reduces to 
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( 43 ) 



and equation (41) can then be expressed as 

{ejf')} = h[A(t)Y'{5 x } = [PWp,}, (44) 

where the time dependence of the laminate is now explicitly indicated. Here we take the applied 

loads to be constant, calculating the time-dependent strain response arising from the time 

dependence of each lamina as characterized by the S-bar (and thus Q-bar) modulus matrices. As 

now defined, the P-bar matrix represents the effective modulus matrix of the laminate in the 

loading coordinate system. For example, using the common notation to let S x denote the 

effective compliance in the x-direction of the laminate 

S x (t) = P u (t)' (45) 

In this manner, the effects of physical aging on the time dependent response of a laminate 
to constant stress states can be determined, using the basic concepts of physical aging to express 
first the transverse and shear compliance of a lamina in the fiber directions, then transforming 
those expressions to obtain the compliance of each lamina in the x-direction, then building up 
the overall laminate response through lamination theory on the compliance of each lamina. As an 
illustration of the results one can obtain with this method, consider a quasi-isotropic laminate, 
[0/±45/90] s . In Figure 16 the change in compliance of this laminate (using the IM7-8320 
composite data) is given over a long time range, plotted on a semi-log scale. Note that even 
though the quasi-isotropic laminate is a fiber dominated lay-up, the compliance changes by 8%- 
12% over a 10 year period. Such a magnitude of modulus change is certainly significant for 
structural applications and must be considered in the design stages. It should be pointed out here 
that the effects of physical aging in terms of stiffness of the composite under long-term loading 
are beneficial: that is, if increase of relaxation times of the matrix was not occurring over time, the 
long time compliance for the composite would collapse to the momentary curve value, also 
indicated for S x on Figure 16. In this case, the compliance increases (and the stiffness decreases) 
at a very rapid rate under constant loading, undesirable for load bearing structures. In view of the 
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section on Aging to Equilibrium, where the momentary response is obtained as a long-term 
response, it is thus also important for composite laminate systems to assure that the equilibrium 
a 8* n g hme will never be reached at the times and temperatures of use. The possible adverse 
effects of physical aging include increased brittleness or decreased fracture toughness with aging 
(Chen-Chi, Lee et al., 1992). 

Conclusion 

The work presented here has outlined a complete analysis of the modeling of physical 
‘tgttig of polymeric composite laminates. The theory for handling aging of polymers was given a 
much needed review and presented in detail. In addition, the theory for aging to equilibrium 
(applicable to polymers and polymer composites alike) was presented for the first time, and a 
new formalization of aging-time and time-temperature relationships was presented. The physical 
aging theory was built up consistently from the pure polymer case, extending to response of a 
composite lamina and then finally a general laminate. The methodology can be applied to easily 
and accurately assess the long-term response of composite systems, considering both the effects 
of different matrix material properties and different composite lay-ups. A critical feature of the 
modeling is that the long-term response is calculated based upon several clearly defined 
parameters that are measured in short-time laboratory tests. 

It was demonstrated that the aging shift rate parameter is the most critical parameter in 
determining magnitude of agings’ effect on long-term response. It was also illustrated in the 
calculation of the response of just a single off-axis lamina, that the interplay of the various time- 
dependent parameters and transformations involved prohibit intuitive guessing of, for example, 
the “least compliant” fiber orientation. This underscores the need for models such as this one to 
be able to reasonably predict the long-term response of composites without extensive long-term 
laboratory testing. It was also demonstrated that even for fiber-dominated lay-ups, physical aging 
is extraordinarily important in calculating the long-term response of laminates accurately. 
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creep compliance (Pa’ 1 ) 



Figure 1 : Evolution of volume towards equilibrium with time below the glass transition 

temperature 



Figure 2a: Momentary creep curves for epoxy with progressing aging time. Experimental 
data to approximately 10 10 Pa 1 for each aging time (Lee and McKenna, 1990). 
Aging to equilibrium phenomenon also illustrated for longer aging times. 
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Figure 2b: Isochronal values of compliance from momentary creep curves in Figure 2a 
versus increasing initial aging times. 



Figure 3: Definition of Shift Rate 




Compliance 



Figure 4: Comparison of long-term and momentary response for polymer with idealized 
properties 
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Creep Compliance 




Figure 6: Bilinear shift rate with equilibrium aging time, L 




t ere f=1680, S 0 =2x10 10 , m=0.276, t£=4.8x10 4 
1 MPa, T=33.2°C, pi=0.858, p 2 =0.12 


t"=1680 

momentary 


t"=1680 
long term 


t°= 1.2x10 s 
momentary 
— long term 


time (sec) 

Momentary and long term response when aged to equilibrium; data after Lee and 
McKenna (1990). 






Figure 8: Schematics of compliance curves at two temperatures and two aging times to 
visualize aging-time and temperature relationship. 



Figure 9: Typical coordinate directions for composite lamina. 
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Compliance in Load direction 
|i22=0.77, p.66=0.93 
T=170°C 



Figure 10: Comparison of the effect of aging on the long term compliance for varying fiber 
angles. Data on IM7/8320 (Hastie and Morris, 1993). 



Figure 1 1 : Long-term compliance for various fiber angles all shifted vertically to 0° curve. 
Illustrates differing effects of aging. 
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Compliance in Load direction 
1122 = 1166 = 0.77 
T=170°C 



Figure 12: Comparison of the effect of aging on the long term compliance for varying fiber 
angles. Data on IM7/8320 (Hastie and Morris, 1993), shift rate modified. 



Figure 1 3: Long-term compliance for various fiber angles all shifted vertically to 0° cun/e. 
Data on IM7/8320 (Hastie and Morris, 1993), shift rate modified. 
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